Post-Newtonian constraints on f(R) cosmologies in Palatini formalism 
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We compute the complete post-Newtonian limit of the Palatini form of f(R) gravities using a 
scalar-tensor representation. By comparing the predictions of these theories with laboratory and 
solar system experiments, we find a set of inequalities that any lagrangian f(R) must satisfy. The 
constraints imposed by those inequalities allow us to find explicit functions that bound from above 
and from below the possible nonlinear terms of the lagrangian. We conclude that the lagrangian 
f(R) must be almost linear in R and that nonlinear corrections that grow at low curvatures are 
incompatible with observations. This result shows that modifications of gravity at very low cosmic 
densities cannot be responsible for the observed cosmic speed-up. 
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I. INTRODUCTION 



The physical mechanism that drives the recently 
observed accelerated expansion rate of the universe 0, 
is yet unclear. Among other possibilities, it has been 
suggested that it could be due to modifications of the 
gravitational interaction at very low cosmic curvatures 
Pi 0] (see also ||). Since the addition of positive 
powers of the scalar curvature to the Hilbert-Einstein 
lagrangian may lead to early-time inflation, it is tempt- 
ing to introduce corrections that grow at low curvatures 
to see if the resulting cosmological models can in some 
way justify the current observations. The theories 
constructed in this way and, with more generality, all 
theories where the gravity lagrangian is an arbitrary 
function of the scalar curvature, are known as f(R) 
gravities. They were originally formulated within the 
standard metric formalism and then exported to the 
Palatini formalism, where metric and connection are 
seen as independent fields. In this work we will study 
the Palatini formulation of such theories. The metric 
formalism will be considered elsewhere la. 



It has been shown that in the Palatini formalism 
many different f(R) lagrangians may lead to late-time 
cosmic acceleration 0, 0, HI- This is true even in 
cases in which the function f(R) is completely different 
from the Hilbcrt-Einstcin lagrangian, i.e., when f(R) 
is not just a perturbation of the linear lagrangian. 
Though much work has been carried out recently to 
understand the cosmological aspects of these theories, 
very little is known about their properties in other 
regimes. The solar system, where the observational 
data are extensive and very precise, represents a more 
suitable and cleaner scenario for testing different aspects 
of the gravitational interaction than the cosmological 
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regime, where much "dirty" astrophysics and uncertain 
data may be involved. In fact, if in addition to a 
modified gravitational dynamics, sources of dark energy 
were acting in the cosmic expansion, it would be very 
difficult to distinguish one effect from the other. In the 
solar system, on the contrary, the matter sources are 
well known and, therefore, the gravitational dynamics 
manifests in its own, with negligible "dark" interferences. 



This work is thus aimed to analyze in detail the 
predictions of the Palatini form of f(R) gravities in 
the weak-field, slow-motion or post-Newtonian regime. 
The Newtonian limit of these theories has been stud- 
ied previously in for the particular lagrangian 
f(R) = R + ujqR 2 , and more recently in and (TJ, 
where a general f(R) lagrangian was considered. Our 
work will try to go a step further obtaining the complete 
post-Newtonian metric and placing explicit constraints 
on the form of the lagrangian. In the computation of the 
post-Newtonian metric, we will use a scalar-tensor repre- 
sentation of these theories. Rather than an unnecessary 
complication or a mathematical trick, the scalar-tensor 
form will be very useful in the interpretation of the 
equations of motion and the identification of the matter 
and the geometrical terms. This identification is neces- 
sary in order to carry out correctly the post-Newtonian 
expansion, as we will point out at due time. 



The paper is organized as follows. We first derive the 
equations of motion and show how to obtain the scalar- 
tensor representation out of them. Then we point out 
some details about the post-Newtonian expansion, com- 
pute the lowest-order corrections of the metric and dis- 
cuss the results (the complete post-Newtonian metric is 
given in the Appendix). Eventually, we will study the 
constraints on the lagrangian f(R) that follow from labo- 
ratory and solar system observational data and will com- 
ment on their consequences in the cosmic regime. We 
end the paper with a brief summary and conclusions. 
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II. EQUATIONS OF MOTION 

The action that defines f(R) gravities in the Palatini 
formalism is the following 



J ^ x V-9f{ R )+ S m\9nu,tpm] (1) 



where S m [g^ V i ipm] represents the matter action, which 
depends on the metric g^ and the matter fields tp rn . The 
scalar R is defined as the contraction R = g^ v R^^T), 
where R^ V (T) has the form of the Ricci tensor 



iW r ) 



fJ.v pX 



(2) 



and rjg 7 is an affine connection independent of g^ v . Vary- 
ing eq.Q) with respect to the metric we obtain 



f'(R)R tlu (T) - -f(R)g, v = n 2 T^ u 



(3) 



where f'(R) = df /dR. Note that the trace of this equa- 
tion 



f'(R)R - 2/(JZ) = k 2 T 



(4) 



implies an algebraic relation between R and T more 
involved than that existing in General Relativity (GR), 
say, R = k 2 T. We will denote by R(T) the algebraic 
solution to eq.|@}- 

The variation of the action with respect to the connec- 
tion gives 



~g{5 x J'g M -\€f'g Xl -\5lf'g px 



o 



(5) 

where /' = f'(R(T)) is also an algebraic function of the 
matter terms. Using an auxiliary tensor t^ v = f'g^, 
eq. (0 can be readily solved |l2( ■ The solution states the 
compatibility between the connection and the metric 
t^ u . In other words, can be written as the Levi-Civita 
connection of 

j-a A 

r /3 7 = — (dptxy + dytxfj - dxtpy) (6) 

Inserting this solution for T^ 7 , written in terms of g^ u 
and /', in eq.© we obtain 



R^{g) - ^g^ R (g) 



>e_ R(T)f> - f 

1 



2(/') 2 
j[V^ v f -g^Df] 



make our notation clearer, since tp V and <7 M „ are con- 
formally related, it follows that R(T) = g^R^T) and 
R{g) — g^R^(g) are related by 

R(T) = R(g) + ^pdxf'd x f jOf (8) 

where, recall, /' = f'(R(T)) is a function of T. 
Introducing the following definitions 



where Rp V {g) and R(g) are computed in terms of the 
Levi-Civita connection of the metric g M „, i.e., they rep- 
resent the usual Ricci tensor and scalar curvature. To 



= /' 

v{<t>) = RWf-m 



(9) 
(10) 



where we have algebraically inverted f'(R) to obtain R = 
R(f')> we can write eq-Q in the scalar-tensor form 

^ij,v\9) — , 1 iiv 2S v 

- ^ (d^d v( j>-\g^4>A + (11) 



1 



(V„V„<£ - g^Ucj)) 



The scalar field <fi satisfies eq.(@J, which in the new nota- 
tion turns into 



2V - 4>V = k 2 T 



(12) 



where V' = dV/d<f>. This representation of the theory 
turns out to be equivalent to an u = —3/2 Brans- 
Dicke-like scalar-tensor theory. For more details on this 
scalar-tensor representation and a different derivation of 
this result see 



epr( 

m 



(7) 



According to ea. l|12|) . the scalar field can be alge- 
braically solved as <fi = </>(T), i.e., it is a non-dynamical 
field. Thus, the effect of the nonlinear lagrangian in the 
Palatini formalism is rather different from its effect in 
the metric formalism 6] . In the metric case, the compat- 
ibility between metric and connection gives rise to ad- 
ditional degrees of freedom in the theory, which can be 
encoded in a dynamical scalar field. In the Palatini for- 
malism, however, the independent connection keeps the 
order of the equations of motion and modifies the way 
matter generates space-time curvature. In other words, 
the right hand side of ea. (|ll|l represents a generalized 
energy-momentum tensor of matter in which the trace T 
plays an enhanced role by means of the terms <f> = </>(T) 
and its derivatives. 



III. POST-NEWTONIAN METRIC 

The complete post-Newtonian limit needs the different 
components of the metric evaluated to the following 
orders 5o0 ~ 0(2) + 0(4), g 0j ~ 0(3), 9lJ ~ 0(2) (see 
Hil)- F° r convenience, we will discuss here only the 
lowest-order corrections, goo ~ 0(2), g^ ~ 0(2), which 
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will be enough to place important constraints on the 
gravity lagrangian. The details of the calculations can be 
found in the Appendix. In our calculations, we will use 
coordinates in which the outer regions of the local system 
are in free fall with respect to the surrounding cosmo- 
logical model. This guarantees that the local metric can 
be made Minkowskian at a given distance far from the 
system and fixes all the boundary conditions necessary 
for this problem. Thus, we will compute the metric as 
a perturbation around the Minkowski background, i.e., 
9^v ~ "Hav + fyui/ ■ Since the scalar field is non-dynamical, 
i.e., it is determined by the local matter distribution, 
no boundary conditions need to be set for it. On the 
other hand, once the solution = 0(T) is obtained, it 
could be expanded to different orders of approximation 
in the post-Newtonian expansion using the fact that for 
a perfect fluid T = —p(l + II — 3P/p) « -p + pO(v 2 /c 2 ), 
where p is the rest-mass density, II is the specific 
energy density (ratio of energy density to rest-mass 
energy), and P is the pressure (see chapter 4 of [T^l. 
In this way one would obtain an expansion of the form 
<f)(T) « <j>(-p) + d T (j}\- P 0{v 2 /c 2 ) + ... However, this 
is an unnecessary complication of our notation and, 
therefore, we will keep <fi(T) exact in our calculations. 
Note that this expansion in post-Newtonian orders 
is different from an expansion around the vacuum 
4>(T) ~ 0(0) + <9W>(0)T + . . . such as the one apparently 
considered in 1(1 and usm g the original f(R) 

representation. In their calculations they expanded 
the function f{R) around a de Sitter background 
characterized by a constant curvature i?o- The fact that 
in the Palatini approach R = R(T) (do not confuse 
R{T) with R(g)), implies that R = R(T = 0). Thus, 
an expansion of f(R) around i?o actually represents 
an expansion around T w 0, which is an expansion 
around the vacuum, not an expansion in post-Newtonian 
orders. Since the functional dependence of f(R) with 
T is a priori unknown, there is no guarantee that such 
expansion around the vacuum can be valid in the range 
from T — up to the typical densities inside planets, 
stars or laboratory-size bodies. In fact, the weak- field 
slow-motion limit does not require low densities but 
not too high matter concentrations and low matter 
velocities, v 2 /c 2 <C 1. Thus, the conclusions reg arding 
the Newtonian limit obtained in [l(| and [111 could 
not be valid. This point will be clarified below in 
detail. In our description in terms of a scalar field, 
4>{T) = f'[R(T)}, the role of the 0(T) terms is clear from 
the very beginning: they represent new contributions 
of the matter sources to the equations of motion. In 
consequence, they must be treated as matter terms and 
expanded in post-Newtonian orders, not around the 



A. Second-order corrections 

For convenience, we introduce a dimensionless quantity 
= 0/00: where 0o = 0(0) is the vacuum reference value, 
and define £1(T) = log(0). In order to get hij diagonal 
and to respect the perturbative description, we find that 
fi(T) must be seen, at least, of order 0(v 2 ). We will 
indicate with a superindex the order of approximation 
of each quantity when necessary. To second order, the 
metric satisfies the equations 



,(2) 
'00 
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K p - 



v(4>) 



20 



> 2 



2 ) 



K 2 p - 



V(0) 



20 



(13) 



(14) 



where we have used the gauge condition 1 hpl — ^h^ k 
dkQ. These equations admit the following solutions 



hod (A x ) 



47T0O 
,.2 



fe) [^W +fl(2) 



47T0O 



d 3 x 



4>\x — x'\ 
, [p( t,x') + V^)/ K 2 
4>\x — x'\ 



In these equations, the local term = log[0(— p)] 

represents a new effect that is not present in the general 
Brans-Dicke-like case u ^ —3/2 (see 6]). The contribu- 
tion due to £1 is identically zero only if f(R) is linear 
(GR and GR plus cosmological constant) and, therefore, 
its presence would imply the nonlinearity of the gravity 
lagrangian. It is worth noting that rather than an 
integrated quantity (cumulative effect), it is directly 
related to the local matter density. In consequence, an 
isolated body will contribute to the exterior space-time 
metric by means of the integral terms of eas. (|15fl and 
(|16f) only. If we now put an object in orbit around the 
first one, the metric at the position of this new body 
will be modified by the local term f2 and by the self- 
gravity of the body. For the moment, we will concentrate 
on the integral terms of eas. (|15fl and i|16|) (isolated body). 

Assuming that the main contribution to the metric in 
the solar system is due to the sun, we can express eqs.JTSJ 
and l|16fl outside the sun as follows 



2G 
2 7 G 



© 



£_ 

r 0o 6 

r 0o 6 



(17) 
(18) 



1 This condition was already used in and also in ll ll . In their 
notation, — dk^i = b^. 
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In these expressions, G and 7 are denned as 

87T0O V M Q. 

M Q - My 



G = 



7 



Mq+M v 



(19) 
(20) 



where Mq = J d 3 x'p(t,x')/(f>, My = k~ 2 f d 3 x'[V - 
V(0)/<^] and Vo = V{(f>o). Since the cosmological con- 
stant term Vq /4>o must be negligible in solar system scales 
in order not to affect the local dynamics, we find a con- 
straint on the function f(R). We need to note that the 
value <po is solution of ea. (|12f> with T = 0. Using that 
equation and the definition of V(<f>) in terms of f(R), it 
follows that Vq = f(Ro), where Rq is solution of eq.(@J 
outside the sun, i.e., Rq = R(T = 0). From these consid- 
erations it follows that 



f(Ro) 



f'(Ro) 



Li « 1 



(21) 



where Ll represents a (Large) length scale the same 
order or greater than the solar system and Rq presum- 
ably is of order the cosmological constant A ~ 10 _53 m~ 2 . 

Let us consider now the observational constraints on G 
and 7. It is well known that in dynamical scalar-tensor 
theories the effective G and 7 depend on two cosmic 
parameters, say, the state of the field, </>q, and the range 
m~ l of its interaction, which are the same for all bodies 
0, 0, 0| . In the non-dynamical situation discussed 
here, G and 7 are not universal quantities, i.e., they are 
not the same for all bodies. According to the definitions 
given above after eq. (J^UJl, two bodies with the same 
Mq do not necessarily have the same value My and, 
therefore, may lead to different values of G and 7. 
This is due to the fact that Mq and My are defined 
as integrals over quantities related to <j> = </>(T), whose 
values depend on the structure and composition of the 
body. Obviously, the experimental evidence suppor ting 
the universality of G and the measurements of 7 ss 1 [lq 
indicate that \My/M@\ <C 1. The only cases in which 
My — correspond to GR and GR plus cosmological 
constant, i.e., those cases in which the lagrangian f(R) is 
linear, or V = Vq =constant. All nonlinear lagrangians 
predict a non trivial potential V{<f) and, therefore, a 
non vanishing My, which may give rise to the effects 
discussed above. Unfortunately, the fact that My is 
given as an integrated quantity does not allow us to 
place any explicit constraint on the form of the function 
f(R)- On the other hand, it is quite disturbing the fact 
that a body with Newtonian mass M N = fd 3 x'p(t,x') 
may yield different values of Mq, G and 7 depending 
on its internal properties. Stated another way, a given 
amount of Newtonian mass may lead to gravitational 
fields of different strengths and dynamical properties. 
Since, as far as we know, effects of this type have not 
been observed in laboratory, we expect a very weak 
dependence of Mq on ip. This is equivalent to saying 



that <f) cannot change too much with the density. 

It is worth noting that with the definitions given above 
for G and 7 and neglecting the cosmological constant 



term, we can write for an isolated body h, 



(2) 
00 



2*7, 



where U represents the Newtonian potential. It is thus 
easy to see that the term (h^) 2 /2 of the complete post- 
Newtonian limit (see the Appendix) leads to the PPN 
parameter 0=1, like in GR. The remaining higher-order 
terms of the metric are all affected by <p. 



B. 



f2' 2 ' contribution 



We will now analyze the effect of the term that 
we omitted above in the case of an isolated massive body. 
As we pointed out, this term must be taken into account 
when a test body is placed within the gravitational field 
of another body. Thus, it must be present in any physical 
situation. Neglecting the cosmological constant contribu- 
tion for simplicity, we can write the metric as follows 



,(2), 

'oo 

,(2), 



h\f{t,x) = [2 7 C/(r) 



(22) 
(23) 



where U(r) = GMq/t is the Newtonian potential gener- 
ated by the massive body and Sl( 2 )(T) = log[(j)(T) / fa] is 
a local term that depends on the matter density T = — p 
at the point (t,x), where our test body is located. First 
thing we need to note is that fi' 2 '(T) is a perturbative 
quantity of order 0(v 2 /c 2 ) <C 1. The only manner to 
respect the perturbative approach is accepting that <fi(T) 
depends very weakly on p, i.e., that 4>(T) must be almost 
constant over a wide range of densities and can be well 
approximated by <f>(T) = O + (90/9T)| T=o T + . . ., 
with 0o 1 (90/9T)| T=o 7 1 < 1 from T = up to nuclear 
densities (T — 10 14 g/cm 3 ) at least. The need for this 
expansion about <po indicates that the lagrangian must 
be almost linear in R (recall that (j> = df/dR). Further- 
more, if <fi(T) had a stronger dependence on T, individual 
atoms could experience strong accelerations due to sud- 
den changes in fl^ when going from outside atoms to 
inside atoms. Those individual microscopic gravitational 
effects would manifest in the macroscopic, averaged, 
description of matter. Since such effects have not been 
observed, they must be very small, if they actually exist. 
Thus, the weak dependence of on T within this wide 
density interval confirms that the contribution of the 
nonlinear terms to the lagrangian f(R) must be very 
small. This conclusion agrees with our previous claims 
regarding the weak dependence of Mq, G and 7 on (f>(T). 

Let us analyze in detail the dependence of <j) on T. We 
will consider, as an illustration, the Newtonian limit of 
the conservation equations V = of a perfect fluid. 
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These equations lead to 



(24) 
(25) 



where a modification with respect to the classical Eulcr 

(2) 

equations is introduced by the term fl contained in h 0Q . 
This modification is given by 



2 



and requires that the condition 

p{d(t>/dT) 



< 1 



(26) 



(27) 



be satisfied over the wide range of densities mentioned 
above in order to guarantee the validity of the macro- 
scopic classical Euler equations. Note that eo.l(2*7)l must 
be true in general, since the contribution of to the ac- 
celeration of a body is given in terms of V^ 2 >. This con- 
straint can be rewritten using eq. (Q2J to evaluate d<t>/dT 
as follows 



(0V" - V) 



< l 



(28) 



It is remarkable the fact that the denominator \<j>V" — V] 
in ea. (|28|l is the counterpart of the effective square mass 
to 2 = [0n Vq" — Vq'] / [3 + 2u>] associated to dynamical 
Brans-Dicke-like fields with u> ^ —3/2 6]. For our discus- 
sion it will be more convenient to see this effective mass 
as an inverse length defining the interaction range of the 
scalar field. We can thus interpret ea. l|28[l as the quotient 
of two lengthscales, one associated to the scalar field over 
another related to the mass density, L~ 2 (p) = (n 2 pc/(f>o). 
Eq. (I28|l can then be seen as the counterpart of the condi- 
tion m^L 2 1 that corresponds to uj ^ —3/2 theories. 
Written in terms of the lagrangian f(R), ea. (|28J) turns 
into 



Rf'(R) 



f'(R) 



Rf"(R) 



L 2 (p) » 1 



(29) 



where /' = /'//g = 4>/4>o- According to our interpre- 
tation of the denominator of eq. l|28|l . we must demand 
that the multiplicative factor in front of L 2 (p) in ca. H29|) 
satisfies 



f'(R) 



Rf"(R) 



> 



l 2 Rf 



(30) 



where I 2 represents a lengthscale much smaller than 
L 2 {p) at nuclear densities, which is the shortest L 2 (p) 
that we can associate to ordinary matter. This inequality 



will allow us to find out which lagrangians satisfy the con- 
dition given in ea. (|29|l at densities such that L 2 (p) ^> I 2 . 
Notice that the modulus of eq. 1)30(1 may lead to more than 
one solution. Let us consider first the case 



f'(R) 



- 1 > 



Rf"(R) ~ l 2 Rf 
Defining A = l/(Rf) > 0, eq.(EH) turns into 



i RdA 

A dR 

1 i RdA 
'A dR 



> 



A 
P 



(31) 



(32) 



Since the right hand side of this equation is positive, the 
left hand side must also be positive. This can only hap- 
pen if the denominator is negative and the numerator is 
positive. In fact, if we denote — e = 1 + ea. (|32|l can 

be written as (1 — e)Je > A/l 2 . Thus, only if e < 1 the 
condition of the sign can be satisfied. In addition, for a 
highly linear lagrangian we expect A/l 2 ^ 1. This leads 
to e <^ 1, which is compatible with The sign and 

magnitude of the denominator indicate that / must be 
very small and positive. A little algebra is enough to 
show that 



d\og[A(A + 2l 2 )} >dlog[i] 



(33) 



Once integrated, the new inequality can be written as 

(A - A + )(A - A_) > (34) 
where an integration constant, Cp, appears in 



A+ = 



~l 2 R±^/c 2 + (l 2 R) 2 
R 



(35) 



Since A must be positive (<f> = /' > to have a well-posed 
theory), the only valid solution to ea. (|34|l is A > A + , 
which implies 



df 

< -f- < 
dR ~ 



(36) 



where we have fixed cq = /q to eliminate the tilde from 
/'. We can finally integrate this last inequality to obtain 



l 2 R 2 R 



fa 



- ,/(/ ) 2 + (l 2 R) 2 +^log[l 2 R+VT+WW] 



21 



f < ot- 

(37) 

Before commenting this result, let us consider the second 
inequality that follows from eq.^SOJ 



f'(R) 



l < 



i 



Rf"{R) ~ l 2 Rf 
Using again the function A = l/(Rf), we obtain 

.4 



9 X RdA 

'A dR. 

1 4- RdA 

A dR. 



> 



I 2 



(38) 



(39) 
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This inequality can only be satisfied if (l + -f^jl) > 0, 
which corresponds to /" < 0. Simple algebraic manipu- 
lations lead to 



df 



> 



dR 

which integrated gives 



f>a- 



l 2 R 2 



(l 2 R) 2 - l 2 R 



(40) 



+YV (/ o )2 + (l 2 R) 2 +^^g[l 2 R+^T+ (I 2 

(41) 

Let us discuss now the significance of eas. H37|) and 
(|41|l . The constant I 2 was introduced in ea. (|30|l to 
represent the length scale over which the nonlinear 
contributions of the gravity lagrangian were relevant. 
For a given I 2 , the nonlinear effects will begin to be 
important about a certain high density scale at which 
L 2 (p)/l 2 ~ 1. If, for instance, we take I 2 = 0, the 
nonlinear effects would be completely suppressed, since 
then ea. (|29[l would be satisfied at all densities. The 
choice I 2 — forces the lagrangian to be linear, which 
can be seen from eas. (|37|l and (|41() in the limit I 2 — ► 0. 
This limit also indicates that the constant /q can 
be naturally set to unity. On the other hand, if the 
nonlinear terms were relevant at low cosmic curvatures, 
the lengthscale I would be of order the radius of the 
universe and, therefore, the nonlinear effects would 
dominate the gravitational dynamics at all scales. 
This fact is obviously in contradiction with our experi- 
ence, as we have discussed in detail throughout this work. 

Though the inequalities derived above are only strictly 
valid in the limit of relatively low curvatures, l 2 R <C 1 
(far from the early-time inflationary period), eas. (|37[l and 
(|4*T|| not only estimate the leading order of the possi- 
ble nonlinear corrections, but give precise functions that 
bound the nonlinearities of the gravity lagrangian in this 
limit. Expanding around l 2 R <C 1 we find 



failure to correctly identify the matter terms and the ge- 
ometrical terms. In any case, expanding O^ 2 -* (T) around 
T = we reproduce the term Ap{x) obtained in [Tlj . 
This term, however, is not present in [T(il ]. On the other 
hand, a Yukawa-type exponential correction in the New- 
tonian potential was found in |lf)j | and |llj |. which is not 
present in our calculations. In the case of a dynami- 
cal field satisfying a second-order differential equation, 
-a term of this type is expected to be related to the in- 
^) 2 ]teraction range of the field. In the Palatini case, the 
field is non-dynamical and, therefore, there is no reason 
for such a term. Moreover, assuming spherical symme- 
try, the Palatini equations admit exact Schwarzschild- 
de Sitter solutions ds 2 = -A(r)dt 2 + dr 2 /A(r) + r 2 dfl 2 
with A(R) = 1 - a/r + Ar 2 /3. The effect of the asymp- 
totic background curvature is given by the Ar 2 /3 term, 
which has the same form as the cosmological constant 
term (Vo/</>o)r 2 that appears in our cas. ((T5|l and IT5) |. 
Thus, there is no reason to expect an exponential cor- 
rection related to the background curvature. The error 
seems to be due to a failure in the identification of the 
leading-order contribution of the term (Rf — f)g^/f- 
In our case, this term is represented by (V/cfyg^v, and 
its leading order is iV/fj))^^ not (V70)/i M „, which could 
justify the Yukawa-type correction for V/(j) =constant. 



l 2 R 2 

R-^<f(R)<a 



R 



l 2 R 2 



(42) 



which confirms that the lagrangian is almost linear in 
R and that the leading order corrections can grow, as 
much, quadratically in R. 



C. On the Newtonian limit 

Before concluding, we will briefly discuss the Newto- 
nian limit obtained in 10] and 11]. As we mentioned 
above, the expansion around the vacuum carried out in 
those papers is not valid a priori. We have shown, how- 
ever, that a viable theory must admit such expansion be- 
cause of the experimental evidence supporting the weak 
dependence of 4> on T. This conclusion, in our case, 
came out after analyzing the predictions of the theory. 
In and [Tlj . the expansion was due to an apparent 



IV. SUMMARY AND CONCLUSIONS 

In this work we have computed the complete post- 
Newtonian metric of the Palatini form of f(R) gravi- 
ties using a scalar-tensor representation of those theories, 
which corresponds to the case ui = —3/2 of Brans-Dicke- 
like scalar-tensor theories. We have also discussed the 
experimental constraints on the lowest-order corrections 
of the resulting metric. We have found that the presence 
of nonlinear terms in the gravity lagrangian makes the 
post-Newtonian metric strongly dependent on the local 
properties of the gravitating system. In particular, the ef- 
fective gravitational mass M , Newton's constant G and 
the PPN parameter 7 of an isolated body may depend on 
the internal structure and composition of the body. This 
follows from the fact that those quantities are given in 
terms of integrals that depend on V (<fi) and <j>, which are 
functions of the trace T. Since those predictions are not 
compatible with observations, we have concluded that 
4>{T) should depend very weakly on T. This point has 
been confirmed by the requirement that the contribution 
f}( 2 ) (x) should be much smaller than unity in order to re- 
spect the perturbative approach. Analyzing the effect of 
the Q( 2 \T) term in the Newtonian limit, we have found 
that the nonlinear corrections at relatively low curva- 
tures (l 2 R <C 1) are bounded by eas. l|57Jl and lj4*T|) . which 
are compatible with an almost linear lagrangian. These 
results show that f(R) gravities in the Palatini formal- 
ism with nonlinear terms that grow at low curvatures, 
such as the Carroll et al. model suggested in Q, 
cannot represent valid mechanisms to justify the cosmic 



7 



acceleration. Nonetheless, Palatini models in which the 
nonlinear terms are suppressed below the limits imposed 
by eqs.(j2IJ) and l|41|l. may still be compatible with ob- 
servations, though they lead to a theory virtually undis- 
tinguishable from General Relativity plus a cosmological 
constant. 
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The contribution coming from the scalar field terms can 
be written as 



(A5) 

where f2 = log(0/0o) and 4>q is an arbitrary constant that 
may be fixed as 0o = 4>(T = 0). The components of t^ v 
are 



2 



'oj 



'00 



n + ^d k nd k h$ 



dodjfl 

2 2 
1 

2 
V 
2d 



(A6) 
(A7) 



K - , - d k n 



h %3 

~2 



-(1 - C) + '-^-dtdjSl + ^h$V 2 Q (A8) 



APPENDIX A: DETAILED CALCULATIONS 

The equations of motion for the metric in Brans-Dickc- 
likc theories are given by 



20 



9fiu p0 + V(0)] 



(Al) 



The expansion of the different components of R^ v can 
be found in the Appendix of 0. We will just remark 
here that in the case oj = —3/2 the scalar field cannot 
be expanded in the same manner as in the general case 
oj =/= —3/2, since now the field is non-dynamical and is 
completely determined by the matter distribution of the 
local system, <f> = <j)(T). In the general case, however, 
the field is a dynamical entity whose state is determined 
by the Universe as a whole. The post-Newtonian system 
only contributes with local fluctuations from the back- 
ground asymptotic state. In the oj = —3/2 case, due to 
the fact that T = -p(l + II - 3P/p) w -p + pO(v 2 ), we 
could expand <j){T) « <j>{—p) + dT<t)pO(v 2 ), though this 
seems an unnecessary complication. We will keep all (j> 
terms exact in our calculations and will expand them at 
the end up to the necessary order. According to this, 
y^T] is given by 



K 2 p 



T 0j 



too 



^S ij+ pO(v 2 ) 
-—pvj + pO(v 3 ) 



2 

up 
~2&~ 



i + n + 2v 2 - 4o + — 



3_P 

p 



(A2) 
(A3) 
pO(w 4 $A4) 



where we assume Q ~ 0(v 2 ) at least to guarantee hij 
diagonal and a consistent post-Newtonian expansion. 

Equating the left hand side of eq. (|A1|) to its right hand 
side, a little algebra leads to 
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(All) 
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Using the gauge conditions 
1 



hit. 



the equations of above become 
1 



(A12) 
(A13) 



2FJ,(2) 

ij ^ «y 



Sufi] 



2 V n oj 4 "-oo,oi 



(A15) 



K 



"-00 " 



^[i + n + 2 , 2 



- ^(1+^3) (A16) 

(2) 

where fyy] denotes the ij-component of h\j and we have 

(2) 

used that, to second order, /i 00 satisfies 



ft (2) -ft 
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by (we denote 
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(A18) 
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(A20) 
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(A21) 



post-Newtonian corrections to the metric are thus given 
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